gusucode.com > 现代通信系统——使用MATLAB(+全部程序) > 现代通信系统——使用MATLAB(+全部程序)/现代通信系统——使用MATLAB(+全部程序)/Matlab 程序/Chapter6/ip_06_08.m
% MATLAB script for Illustrative Problem 6.8. echo on N=31; T=1; W=1/(2*T); n=-(N-1)/2:(N-1)/2; % the indices for g_T % The expression for g_T is obtained next. for i=1:length(n), g_T(i)=0; for m=-(N-1)/2:(N-1)/2, if ( abs((4*m)/(N*T)) <= W ), g_T(i)=g_T(i)+sqrt((1/W)*cos((2*pi*m)/(N*T*W)))*exp(j*2*pi*m*n(i)/N); end; echo off ; end; end; echo on ; g_T=real(g_T) ; % The imaginary part is due to the finite machine precision % Obtain g_T(n-(N-1)/2). n2=0:N-1; % Obtain the frequency response characteristics. [G_T,W]=freqz(g_T,1); % normalized magnitude response magG_T_in_dB=20*log10(abs(G_T)/max(abs(G_T))); % impulse response of the cascade of the transmitter and the receiver filters g_R=g_T; imp_resp_of_cascade=conv(g_R,g_T); % Plotting commands follow.